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ABSTRACT 

We present the analysis of the temperature and metallicity profiles of 12 galaxy clusters in the 
redshift range 0.1-0.3 selected from the Chandra archive with at least ~ 20, 000 net ACIS counts and 
kT > 6 keV. We divide the sample between 7 Cooling-Core (CC) and 5 Non- Cooling- Core (NCC) 
clusters according to their central cooling time. We find that single power-laws can describe properly 
both the temperature and metallicity profiles at radii larger than O.lrigo in both CC and NCC systems, 
showing the NCC objects steeper profiles outwards. A significant deviation is only present in the inner 
O.lrigo. We perform a comparison of our sample with the De Grandi & Molendi BeppoSAX sample of 
local CC and NCC clusters, finding a complete agreement in the CC cluster profile and a marginally 
higher value (at ~ la) in the inner regions of the NCC clusters. The slope of the power-law describing 
kT{r) within O.lrigo correlates strongly with the ratio between the cooling time and the age of the 
Universe at the cluster redshift, being the slope > and T c /T age < 0.6 in CC systems. 

Subject headings: galaxies: clusters: general — (galaxies:) intergalactic medium — X-rays: galaxies: 
clusters 



1. INTRODUCTION 

Clusters of galaxies represent unique signposts in the 
Universe, where the physical properties of the cosmic dif- 
fuse baryons can be studied in great details and used 
to trace the past history of cosmic structure formation 
(e.g. Rosati et al. 2002; Voit 2005, for reviews). As 
a result of adiabatic compression and shocks generated 
by supersonic motion during shell crossing and virializa- 
tion, a hot thin gas permeating the cluster gravitational 
potential well is formed. Typically this gas, which is en- 
riched with metals ejected form supernovae (SNe) explo- 
sions through subsequent episodes of star formation (e.g. 
Matteucci & Vettolani 1988; Renzini 1997), reaches tem- 
peratures of several 10 7 K and therefore emits mainly via 
thermal bremsstrahlung in the X-rays. At such temper- 
atures most of the elements are either fully ionized or in 
a high ionization state. 

Particularly evident in X-ray spectra of galaxy clusters 
are the strong transitions to the n = 1 level (K-shell) of 
the H-like and He-like ions of Iron in the energy range 
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6.7-6.9 keV. Below 2 keV, n = 2 level (L-shell) transi- 
tion of Iron and a elements can be detected, especially in 
the low temperature region in the centers of the so-called 
cool core clusters, that are characterized by a strong peak 
in the surface brightness distribution, and therefore short 
cooling times. Spatially resolved CC clusters show a peak 
in the metal distribution associated to the low temper- 
ature core region (e.g. De Grandi & Molendi 2001 and 
2002, hereafter DM01 and DM02). Although the amount 
of energy supplied to the intra-cluster medium (ICM) 
by SNe explosions depends on several factors (e.g. the 
physical condition of the ICM at the epoch of the en- 
richment) and cannot be obtained directly from X-ray 
observations, the radial distribution of metals, as well as 
their abundance as a function of time, are crucial infor- 
mation to shed light on the cosmic star formation history 
and to trace the effect of SN feedback on the ICM. 
Several analyses have been presented in the literature 
with the aim to study the radial distribution of met- 
als in clusters of galaxies. Finoguenov, David & Ponman 
(2000) performed a spatially resolved X-ray spectroscopic 
analysis of 11 relaxed clusters observed by ROSAT and 
ASCA, deriving a radial distribution of single heavy ele- 
ments such as Fe, Si, Ne and S. They found that the to- 
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tal Fe abundance decreases significantly with radius in all 
clusters, while the Si, Ne and S abundances are either flat 
or decrease less rapidly. DM0 f derived radial metallicity 
profiles (mainly driven by Fe) of 17 nearby clusters ob- 
served by BeppoSAX. They found a strong enhancement 
in the abundance in the central regions of the CC clus- 
ters. A flatter metallicity profile was observed instead for 
the non-cool core clusters in their sample. Since all the 
NCC clusters show signs of recent merger activity, they 
suggested that the merger events may have redistributed 
efficiently the metal content of the intracluster medium. 
Irwin & Bregman (2001) derived iron-abundance pro- 
files for 12 clusters with 0.03 < z < 0.2 observed by 
Beppo SAX. Although they investigated the differences 
between CC and NCC clusters in a less systematic way 
than DM01, they found a negative gradient in the abun- 
dance profiles of all the CC clusters and to a lesser sig- 
nificance also in the NCC clusters. Similarly to DM01, 
they found that CC clusters have higher metallicity than 
NCC clusters at every radius. It is worth to say that 
the aforementioned papers investigated the metallicity 
trends only within rsoo- Spatially resolved measures of 
the metal abundance in galaxy clusters were performed 
also with XMM-Newton. In particular Tamura et al. 
(2004) analyzed a sample of 19 X-ray bright relaxed clus- 
ters, obtaining elemental abundances of Fe, Si, S and O. 
They found that while the distribution of Fe, Si and S 
is generally peaked toward the center, the O abundances 
are uniform throughout the cluster, pointing out to a dif- 
ferent origin among these metals, most likely in SNe la 
and II. More recently Vikhlinin et al. (2005) have derived 
temperature and metallicity profiles for 11 low-redshift 
clusters observed by Chandra. The clusters in their sam- 
ple however are all CC clusters, presenting a very regular 
overall X-ray morphology and showing only very weak 
signs of dynamical activity. Although they have not an- 
alyzed the metallicity profiles rescaled to the virial radius 
of the cluster (as they did for the temperature profiles) a 
negative gradient of Z is present in all the objects in their 
sample. However, almost all the spatially-resolved metal- 
licity profiles are measured in local clusters (z < 0.1). 

On the other hand, measurements of the metal content 
of the intracluster medium at high— z has been obtained 
with single emission-weighted estimates from Chandra 
and XMM-Newton exposures of 56 clusters at 0.3 < z < 
1.3 in Balestra et al. (2007). They measured the Iron 
abundance within (0.15-0. 3) R V i r and found a negative 
evolution of Zp e with the redsfhift, with clusters at z > 
0.5 having a constant average Fe abundance of f» 0.25Z Q , 
while objects in the redshift range 0.3 < z < 0.5 show 
Zpe significantly higher (« OAZq). Such evolution is not 
driven entirely by the presence of the cool cores. This 
result has been recently confirmed by Maughan et al. 
(2007). 

In this paper, we present measurements of the radial 
temperature and metallicity profiles of a sample of 12 
clusters with temperatures larger than 6 keV observed 
with Chandra at intermediate redshift, 0.11 < z < 0.32. 
We take advantage of the ACIS superior spatial and spec- 
tral resolution to investigate in a systematic fashion the 
differences that may exist between CC and NCC clus- 
ters. The spectroscopic measurements of the ICM tem- 
perature and metallicity allow to characterize statisti- 
cally the radial profiles and to quantify their gradients in 



this unexplored redshift region. 

All the uncertainties are quoted at la (68%) for one 
interesting parameter. The abundance estimates are rel- 
ative to the compilation of cosmic values given in Anders 
& Grevesse (1989) (hereafter AG89), unless otherwise 
stated. Indeed, these values for the solar metallicities 
have more recently been superseeded by the new values 
by Grevesse & Sauval (1998) and Asplund et al. (2005) 
(hereafter A05), who introduced a 0.676 and 0.60 times 
lower Iron solar abundance, respectively (photospheric 
value), while the other elements do not change signifi- 
cantly. Our measures of metallicity are expected to be 
driven mainly by Iron, however, for clarity, we also per- 
formed the fits using solar abundances by A05. 
Throughout this paper we assume Hq = 100 h km s _1 
Mpc" 1 , h = 0.7, Q m = 0.3 and n A = 0.7. 

2. SAMPLE DEFINITION AND DATA ANALYSIS 

From Chandra archival data we select a sample of 
twelve 'intermediate' redshift clusters (0.11 < z < 0.32). 
We also require the clusters to have at least ~ 20, 000 
ACIS-S or ACIS-I counts in order to study their 
properties in at least 3 circular annuli. The sample is 
presented in Table [H where the name of the cluster and 
the Chandra observing logs are listed. 

The Chandra data analysis has been performed 
using the latest version of CIAO (v3.3.0.1). All of our 
datasets are processed by a version of the Standard 
Data Processing (SDP) pipeline prior to version DS 
7.4.0, which uses the tool acis-detectjxfterglow to 
flag possible cosmic ray events in the level 1 event 
file; it has been determined that a significant fraction 
of the X-ray events from a source in imaging mode 
might be removed using this tool. Therefore we reset 
the correction performed by acisjdetect_afterglow 
on the Level=l event file, so that the hot pixels and 
the afterglow events may be properly removed by the 
improved CIAO tool acisjrunJiotpix, (introduced after 
SDP version DS 7.4.0). A new Level=l event file is then 
created (through the CIAO tool acis-process-events) to 
apply the latest calibration files to the data (e.g. apply 
the newest ACIS gain maps, apply the time-dependent 
ACIS gain correction, apply the ACIS Charge Transfer 
Inefficiency correction, etc.). Moreover in the case of 
observations telemetered in VFAINT mode it is possible 
to reduce the background using the additional screening 
of the events with significantly positive pixels at the 
border of the 5x5 event island. Two further filtering 
steps are then required to obtain the Level=2 event 
files, i.e. filter for bad grades (using ASCA grades) and 
for a "clean" status column and apply the Good Time 
Intervals (GTIs) supplied by the pipeline. The final 
step is to examine background light curves during each 
observation to detect and remove the periods of high 
background, due to flaring episodes. We perform the 
flare detection and removal following the recommen- 
dations suggested in Markevitch et al. (2003); both 
the total and the clean exposure times are listed in 
Table [TJ Most of the observations are slightly affected 
by background flares thus we were able to use practically 
all the exposure time. The only exceptions are ObsID 
3205, 4193, 906 and especially ObsID 1666, where - 29 
ks of the exposure were lost due to high background. 
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TABLE 1 

Observation Log for the Chandra Cluster Archive Sample. Column (2) shows the redshift z of the clusters. The columns 
(3), (4) and (5) show the instrument used, the observation date and the observation id, respectively. the columns (6) and 

(7) SHOW THE OBSERVING TIME BEFORE (t exp ) AND AFTER (t c i ean ) THE REMOVAL OF HIGH BACKGROUND INTERVALS. COLUMN (9) IS THE 

Galactic column density N h in the line of sight of the observation. 



Name 


z 


Instrument 


Obs. Date 


Obs. ID 


texp 

(ksec) 


^ clean 

(ksec) 


N H 
(10 20 cm" 2 ) 


A2U34 


0.113 


A PTC T 


2UU1 May UO 


2204 


53.9 


53.9 


1.6 


A1413 


0.143 


ACIS-I 


20U1 May 16 


1661 


9.7 


9.7 


2.2 






ACIS-I 


2004 Mar 06 


5003 


76.1 


75.0 








ACIS-I 


2005 Feb 03 


5002 


37.2 


36.5 




A907 


0.153 


ACIS-I 


2000 Jun 29 


535 


11.0 


10.9 


5.4 






ACIS-I 


2002 Jun 14 


3185 


48.7 


47.9 








ACIS-I 


2002 Oct 30 


3205 


47.7 


40.5 




A2104 


0.155 


ACIS-S 


2000 May 25 


895 


49.8 


48.9 


8.7 


A1914 


0.171 


ACIS-I 


2003 Sep 03 


3593 


18.9 


18.8 


0.9 


A2218 


0.176 


ACIS-I 


2001 Aug 30 


1666 


49.2 


20.2 


3.2 


A963 


0.206 


ACIS-S 


2000 Oct 11 


903 


36.8 


35.8 


1.4 


A2261 


0.224 


ACIS-I 


2004 Jan 14 


5007 


24.6 


24.3 


3.3 


A2390 


0.228 


ACIS-S 


2000 Oct 08 


500 


9.8 


9.8 


6.8 






ACIS-S 


2003 Sep 11 


4193 


96.3 


91.0 




A1835 


0.253 


ACIS-S 


2000 Apr 29 


496 


10.8 


10.3 


2.3 


ZwC13146 


0.291 


ACIS-I 


2000 May 10 


909 


46.6 


45.6 


3.0 


A1995 


0.319 


ACIS-S 


2000 May 08 


906 


57.5 


53.8 


1.4 



2.1. Background Subtraction 

An accurate subtraction of the background is crucial to 
perform a correct study of the spectral properties of the 
clusters in our sample, especially in their outskirts. Since 
we are dealing with extended objects, occupying most of 
the ACIS field of view we need to use a compilation of 
the blank-field observations, processed identically to the 
cluster observation (i.e. as described above) and repro- 
jected onto the sky using the aspect information from the 
cluster pointing. It is worth noticing that the synthetic 
backgrounds correspond to longer exposure times (~ 0.5 
Msec) than any of our observations, giving us a very good 
sampling in the estimate of the background to subtract. 
Moreover, in order to 'tailor' the background to our data 
we follow the recommendations given in the CIAO web- 
pages . In particular we renormalize the blank- fields to 
the background in each observation, considering a region 
of the ACIS field of view practically free from cluster 
emission (mainly ACIS-S 1 for ACIS-S observations, and 
ACIS-S2 for ACIS-I observations) and a spectral band 
(9.5-12 keV) where the Chandra effective area is nearly 
zero, therefore all the observed flux is due to the particle 
background. 

In addition to the particle-induced background we 
check also if the diffuse soft X-ray background could be 
an important factor in our observations and if appropri- 
ate adjustments are needed. For each observation, we 
follow the procedure of Vikhlinin et al. (2005), extract- 
ing a spectra in the source- free regions of the detector, 
subtracting the renormalized blank-field background and 
fitting the residuals in XSPEC vll.3.2p (in the 0.4-1 keV 
band) with an unabsorbed mekal model, whose normal- 
ization was allowed to be negative. The best-fit model 
obtained is therefore included as an additional compo- 
nent in the spectral fits (with its normalization scaled by 
the area). However, in every observation the adjustments 

4 http://cxc.harvard.edu/cal/Acis/Cal_proda/bkgrnd/aciabg/COOKBOOK 



required are minimal and do not affect significantly the 
determination of kT and Z, even at large radii. This 
is also due to the properties of the clusters in our sam- 
ple, whose high kT values (> 4 keV) even in the re- 
solved outer regions are not affected significantly from 
the method applied for the subtraction of the diffuse soft 
background. 

2.2. Cash statistics vs. x 2 statistics 

The x 2 statistics require grouping of the spectra, hav- 
ing at least 20 counts per bin, in order to be able to 
approximate the Poissonian distribution of counts with 
a Gaussian. On the contrary, Cash statistics do not re- 
quire any grouping and represent a more reliable (and 
less biased) approach to fit the data. Indeed, it is well 
known from the literature (see e.g. Nousek & Shue 1989; 
Balestra et al. 2007) that the % 2 statistics systemat- 
ically "sees" the observed spectra softer than the real 
ones. This usually leads to an overestimate of the slope 
of the observed spectra in the case of a simple power- 
law fit, while in the case that a thermal model is fitted 
to the data, the temperature measured is usually under- 
estimated. As a test to see whether this systematics is 
present also in our data, we have decided to apply both 
these fit statistics. We find that for all the clusters of 
the sample, a systematically lower temperature is mea- 
sured with the x 2 ( on average 7% — 17% lower, depending 
on the cluster). On the other hand no obvious system- 
atic trend is observed in the determination of Z , being 
the variation in the best-fitting value of the metallicity 
in each cluster ranging between AZ ~ 0.01 and ~ 0.07 
with no preferential direction. To avoid the dependence 
on the grouping method, and the bias in the best-fit tem- 
perature, we decide to use the modified Cash statistics, 
as implemented in XSPEC vll.3.2p, to determine the 
best-fit parameters and their uncertainties. 

2.3. Spectral Analysis 
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In order to study the radial properties of the cluster 
emission, we subdivide each cluster in annuli (circular 
or elliptical, depending on the morphology of the clus- 
ter) centered on the X-ray emission peak. In the more 
disturbed clusters, where an emission peak is not clearly 
identifiable, we assume the center of the cluster to cor- 
respond with the X-ray centroid at 0.5r5oo- We require 
each region to have at least ~ 7, 000 net counts, so that 
it would be possible to estimate the temperature and 
the metallicity of the annulus with sufficient accuracy. 
For each cluster the outermost annulus corresponds to 
an area where the intensity of the source in the 0.8-8 
keV band is roughly equal to that of the background. 
We extract a spectra from each annulus after exclud- 
ing the 3cr point sources detected by the CIAO tool 
wavdetect. The source list produced is also inspected 
'by-eye' in order to remove possible additional sources 
not detected by wavdetect (especially in the regions 
where the diffuse emission from the cluster is brighter). 
The CIAO script used to perform the spectral extraction 
is specextract, which generates source and background 
spectra and build the appropriate RMFs and ARFs. The 
background is taken from the re-normalized blank field 
observations using the same region of the source. 
The spectra are analyzed with XSPEC vll.3.2p (Arnaud 
et al. 1996) and fitted by a single-temperature mekal 
model (Kaastra 1992; Liedahl et al. 1995) in which the 
ratio between the elements is fixed to the solar value as in 
AG89. However, as explained in §1, these values for the 
solar metallicities have more recently been superseeded 
by the new values by Grevesse & Sauval (1998) and A05. 
For clarity and completeness, we also performed the fits 
using solar abundances by A05. The free parameters in 
the model are the temperature kT, the metallicity Z of 
the gas and the normalization. The spectral band consid- 
ered in the fit is the 0.6-8 keV. We choose not to consider 
the data below 0.6 keV because of uncertainties in the 
ACIS calibration below that energy. The Nh derived 
from the X-rays is found to be consistent (within la) 
with the Galactic value in the line of sight of each obser- 
vation, as derived from radio data (Stark et al. 1992), 
except in the cases of A2104 and A2390 (see Tableland 
following § 12. 4|) . In these clusters the Nh value mea- 
sured from X-ray data is significantly different (at more 
than 2cr confidence level) from the radio value, therefore 
we adopt the X-ray value. The Njj value is fixed to the 
Galactic value obtained from the radio data (and listed 
in Table [T]) in the rest of the sample. We have measured 
Nh from the X-ray data in each annulus, finding no ev- 
idence of radial variation. Therefore, Nh is fixed to the 
same value in all the radial annuli. 

We have divided then our sample in Cooling-Core 
and Non-Cooling-Core clusters according to their central 
cooling time. The gas temperature and density profiles 
are recovered from the single-phase spectral fit done in 
annular rings by correcting the emissivity in each shell by 
the contrubution of the outer shells moving inwards. A 
detailed description of the procedure is presented in Et- 
tori et al. (2002). In brief, the normalization of the ther- 
mal component, being proportional to the Emission In- 
tegral, provides the gas density, whereas the deprojected 
temperature is provided by weighting for the corrected 
emissivity the spectral measurement. The deprojected 
values in the innermost bin are then used to estimate 
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Fig. 1. — Histograms of the distribution in our sample of the 
values of T c /T age , computed as described in § [3] and quoted in 
Tabled 



the central cooling times t c = 5/2(/i e //i)T e (n e /e), where 
fi = 0.613 and /i e = 1.174 are appropriate for a plasma 
with a metallicity of 0.3 times the solar values in AG89, 
T e , n e and e are the gas temperature, electron density 
and emissivity in the innermost bin, respectively. 

The central cooling times are reported in Table [2 as 
well as the age r oge of each cluster, and the ratio be- 
tween the two quantities. The age of the Universe at the 
z of observation is used as an upper limit to the age of 
the cluster. Bauer et al. (2005) computed the cooling 
times for 6 of the clusters in our sample (A1835, A1914, 
A2218, A2261, A2390 and ZwC13146) finding a r c in the 
center of the cluster, or at 50 kpc, consistent with the val- 
ues computed for the central bin in our spectral analysis 
(which might extend farther out than 50 kpc from the 
center in some cases). Following their criterium, a clear 
separation between the CC and the NCC in our sample 
can be located at r c ~ 10 Gyr (Figure [TJ correspond- 
ing to r c /r age ^ 1). We have 4 clusters presenting signs 
of strong cooling (r c < 2 Gyr) and 3 clusters exhibiting 
signs of mild cooling (r c < 10 Gyr). The remaining 5 
clusters can be classified as NCC, presenting longer cool- 
ing times in the center. The projected temperature and 
metal abundance profiles for both CC and NCC objects 
are shown in Figures [2] and [3) 

2.4. Notes on individual clusters 

A2034: A2034 (z=0.113) has been observed with 
Chandra in one ACIS-I pointing (ObsID: 2204). The 
temperature profile we derived is quite flat in the central 
regions of the cluster, where the temperature is kT ~ 8 
keV. It shows however a negative gradient after 400 kpc 
from the center. A similar trend is observed in the metal- 
licity, where the average value of Z = 0.4 within 400 kpc 
from the center decreases to Z < 0.2 at larger radii. 

A1413: A1413 (z=0.143) has been observed four times 
with ACIS-I. We discard one observation (ObsID: 537) 
which is affected almost entirely by a persisting flare. 
In one of the observations used in our analysis (ObsID: 
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TABLE 2 

Global cluster properties of the Chandra sample. Column (2) shows the the total net counts from the inner to the 
outer annulus considered in the spectral analysis. the columns (3), (4) and (5) show the global temperature (kt) , the 
global metallicity (z) and the virial radius rigo, computed within 0.07 — 0.4rigo , respectively. column (6) shows the 
aperture used to measure (kt) and (z) . the columns (7) and (8) show the central cooling times and the ratio with 

respect to the age of the universe at the cluster redshift. 



Name 


Net cts 
(0.6-8 keV) 


(kT) 
(keV) 


(Z) 
(Zq) 


riso 
(kpc) 


Aperture 
(") 


To (Gyr) 


7~cf 7~a,ge 


A2034 


78,900 


6.36 ± 0.15 


0.30 ± 0.04 


2,222 


76-433 


21.2 ±3.2 


1.77 ±0.27 


A1413 


181,500 


' - oz -0.12 


0.23 ±0.03 


2,416 


67-385 


4.2 ±0.3 


0.36 ±0.03 


A907 


87,500 


5.82 ± 0.12 


0.34 ±0.04 


2,125 


56-320 


2.0 ±0.1 


0.17 ±0.01 


A2104 


63,100 


6.76 ±0.19 


0.24 ± 0.05 


2,290 


60-341 


18.1 ± 2.1 


1.58 ±0.18 


A1914 


39,100 


9 20+ 039 


0.27 ±0.07 


2,672 


64-367 


12.2 ± 1.0 


1.07 ±0.09 


A2218 


18,300 


6.25 ±0.31 


0.24 ±0.07 


2,202 


52-295 


21.3 ± 1.8 


1.89 ±0.16 


A963 


41,800 


fi 02 +0 ' 28 
D - uz _0.19 


0.18 ±0.06 


2,161 


45-256 


6.5 ±0.4 


0.59 ±0.04 


A2261 


21,500 


7 43+0.49 
' 0.27 


30+ 07 


2,400 


47-267 


7.0 ±0.4 


0.65 ±0.03 


A2390 


202,600 


9.35 ± 0.15 


0.30 ±0.03 


2,693 


52-295 


1.3 ±0.2 


0.12 ±0.02 


A1835 


23,100 


8.06 ±0.53 


0.31 ±0.09 


2,500 


44-254 


0.9 ±0.1 


0.08 ±0.01 


ZwC13146 


40,500 


8.59 ±0.39 


0.24 ±0.06 


2,582 


41-237 


1.0 ±0.1 


0.10 ±0.01 


A1995 


30,200 




0.501°;^ 


2,427 


37-209 


12.7 ± 1.3 


1.28 ±0.14 



5003) the source is placed in a position of the ACIS-I 
array very close to the S2 chip, therefore S2 is still con- 
taminated by source emission and we can not use it to 
re-normalize the blank field to the background in the ob- 
servation. We use instead part of the II chip (which is 
front-illuminated as S2) to re-normalize, since it is more 
distant from the cluster center than S2 and therefore less 
contaminated by cluster emission. The resulting tem- 
perature profile shows a slight decrease in temperature 
towards the center (AkT = 1.2+q'q keV within the in- 
ner 150 kpc). A1413 has been also observed by XMM- 
Newton (Pratt & Arnaud 2002), representing one of the 
clusters with the most accurate temperature profile ob- 
served by this satellite. This cluster is also part of the 
sample of Chandra clusters analyzied by Vikhlinin et al. 
(2005). The XMM-Newton observation does not find any 
evidence of a cool core, in contrast with the temperature 
profiles obtained with Chandra both in our analysis and 
even more evidently in Vikhlinin et al. (2005). This 
might be due to the poorer angular resolution of XMM- 
Newton with respect to Chandra. 

The metallicity profile is decreasing towards larger radii, 
and consistent within la with the measures of Vikhlinin 
et al. (2005). 

A907: A907(z=0.153) has been observed with Chan- 
dra in three separate ACIS-I pointings (ObsID: 535, 3185 
and 3205), all of them used in our analysis. The tem- 
perature profile shows an evidence of a cool core in the 
center of the cluster (AkT — 1.4ftg | keV m the central 
100 kpc). The metallicity profile presents a decreasing 
trend toward larger radii. A907 is also part of the clus- 
ter sample analyzed by Vikhlinin et al. (2005). Their 
results, both for the temperature and the metallicity, are 
fully consistent with ours within the la statistical uncer- 
tainties. 

A2104: A2104 (z=0.155) has been observed with 
Chandra in one ACIS-S pointing (ObsID: 895). As de- 
scribed in § I2.3[ the value of the Nh measured from the 
X-ray data alone is significantly different from the radio 
value, thus we have decided to fix the Nh to the best 
fit value obtained from the fit (1.55 x 10 21 cm -2 ). The 



cluster does not show any evidence of a cool core in its 
center, having a temperature profile decreasing towards 
the outskirts. The metallicity profile is consistent with 
being flat, with a value Z ~ 0.3 — 0.4Zq, within the la 
uncertainties. 

A1914: Two ACIS-I pointings of A1914 (z=0.171) are 
available in Chandra archive. However the oldest (and 
shortest) observation (ObsID: 542) has been performed 
in 1999. For observations performed in that year an accu- 
rate modeling of the ACIS background is not currently 
available. Thus to avoid problems in background sub- 
traction we have decided to discard it and keep only the 
longest observation (ObsID: 3593). A negative gradi- 
ent in kT is quite clear: the temperature drops from 
kT = 12.0;^ keV in the center down to kT = 8.5 ± 0.6 
keV in the outer radial bin. A similar trend is observed 
also in the abundance profile, where Z = 0.5 ± 0.1 Zq 
in the center, then decreasing to 0.2 ± O.IZq in the two 
outer radial bins. This is one of the few examples of 
metallicity peak without a corresponding cool core (or 
temperature drop) towards the center. 

A2218: A2218 (z=Q.176) has been observed three 
times with ACIS-S. Unfortunately two of these observa- 
tions (ObsID: 553 and 1454) were performed in 1999 and 
for the reason described in the case of A1914 we have 
decided to discard them. Moreover the remaining ob- 
servation (ObsID: 1666) has been strongly affected by a 
flare which reduces the good exposure time to only ~ 20 
ks. With these data we are able to observe the presence 
of a centrally peaked temperature profile (a hot, instead 
of a cool core) and a costant metallicity profile. A tem- 
perature profile peaked toward the center has been also 
seen by Machacek et al. (2002), analyzing the two Chan- 
dra observations performed in 1999. This is consistent 
with the picture of A2218 being involved in a line-of- 
sight merger, as suggested by a considerable disturbance 
of the intracluster gas in the X-rays and by the observed 
substructure in the optical (e.g. Pratt et al. 2005). 

A963: A963 (z=0.206) has been observed with Chan- 
dra in one ACIS-S pointing (ObsID: 903). We found a 
decreasing trend of Z with the radius, with only a very 
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Fig. 2. — a): Temperature profiles for the CC clusters in the sample: A907 (green empty triangles), A963 (blue empty squares), Ai413 
(red empty circles), Af835 (magenta filled triangles), A226i (green stars), A2390 (cyan filled circles) and ZwC13146 (red filled squares). 
b): Normalized temperature profiles for the CC clusters, plotted against the radii in units of riso- The symbols have the same meaning as 
in panel a), c): Temperature profiles for the NCC clusters in the sample: Ai9f4 (red filled circles), A1995 (green filled triangles), A2034 
(cyan empty circles), A2f04 (blue filled squares) and A2218 (magenta empty triangles), d): Normalized temperature profiles for the NCC 
clusters, plotted against the radii in units of riso- The symbols have the same meaning as in panel c). 



weak hint for the presence of a lower temperature in the 
center. 

A2261: Two pointings of A2261 (z=0.224) are avail- 
able in the Chandra archive. One of the observations 
(ObsID: 550) has been performed in 1999 and therefore 
we discard it for the reason described above in the case 
of A1914. The temperature profile does show only a hint 
(more than 2a however) of a decrease in the center, where 
the temperature drops down from 9.0±0.4 keV to 7.7±0.4 
keV. The metallicity profile shows a constant behaviour 
for the first two bins and a decrease (significant at more 
than la) in the outer radial bin. 

A2390: A2390 (z=0.228) has been observed three 
times with ACIS-S. One of the observations (ObsID: 501) 
has been performed in 1999 and therefore we discard it. 



We concentrate our analysis on the remaining two obser- 
vations (ObsID: 500 and 4193), yielding a total of ~ 100 
ks of good observing time. The value of the Nh derived 
from the X-rays (Nh = 1.1 X 10 21 cm~ 2 ) is significantly 
different than the radio value, therefore we adopted the 
X-ray value in the spectral fits. A2390 is also part of the 
sample analyzed in Vikhlinin et al. (2005). Similarly to 
them we find a cool core (kT = 5.8 ± 0.2 keV) in the cen- 
ter of the cluster with a kT profile getting flatter going 
towards the outskirts, being fully consistent with their 
measured temperatures at every radius. On the other 
hand, the metallicity profile shows just a hint of a peak 
in the central part of the cluster. It is however consistent 
at la with the profile of Vikhlinin et al. (2005) and not 
sensitive to the choice of the Nh- 
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Fig. 3. — a): Abundance profiles for the CC clusters in the sample: the different symbols correspond to the clusters listed in Figure |2ji. 
b): Normalized abundance profiles for the CC clusters, plotted against the radii in units of riso- The symbols have the same meaning as in 
panel a), c): Abundance profiles for the NCC clusters in the sample: the different symbols correspond to the clusters listed in Figure|2]:. 
d): Normalized abundance profiles for the NCC clusters, plotted against the radii in units of rigo- The symbols have the same meaning as 
in panel c) 



A1835: Two different ACIS-S observations of Af835 
(z=0.253) are present in the Chandra archive. We dis- 
card the older (and longer) observation (ObsID: 495), 
performed in 1999 because of the reason described in the 
case of A1914, keeping only the ~ 10 ks observation per- 
formed in 2000 (ObsID: 496). The temperature profile 
of A 183 5 shows a clear evidence of a cool core in its cen- 
ter where kT drops down by a factor of ~ 2. Moreover 
the temperature shows a decline after 300 kpc, going to- 
ward larger radii. Piffaretti et al. (2005) analysed XMM- 
Newton observations of A1835 and detected a tempera- 
ture decrease at large radii, as in our data. Majerowicz et 
al. (2002) also analysed XMM-Newton data and found 
a decrease in the temperature profile at large radii (at 
~ 400 kpc from the center); however their temperature 
profile becomes constant after such decrease. The de- 



crease at large radii has not been observed in the analysis 
of Chandra data by Voigt & Fabian (2006), who found a 
constant temperature outside the central 100 kpc. How- 
ever it is worth noticing that in their work they analysed 
the 1999 observation (instead of the 2000 observation, as 
in our analysis) which might have background subtrac- 
tion problems especially at large radii. This may explain 
the difference between the two profiles. The metallic- 
ity profile shows a decreasing gradient in the first two 
bins, becoming constant afterwards. The only compar- 
ison with the literature comes from an XMM-Newton 
observation analyzed by Majerowicz et al. (2002), where 
an almost constant metallicity profile at every radius has 
been observed. 

ZwC13146: ZwC13146 (z=0.291) has been observed 
with ACIS-I in one pointing (ObsID: 909). Also this 
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cluster clearly shows the presence of a cool core (kT drop- 
ping down by a factor of almost 2). A decreasing trend 
in metallicity from Z = 0.50 ± 0.05 in the center, down 
to Z = 0.17 ± 0.12 in the outer bin, is also observed. 

A1995: One ACIS-S observation of A1995, the far- 
thest cluster in our sample (z=0.319), is present in the 
Chandra archive (ObsID: 906). Although this observa- 
tion is quite long (~ 50 ks of good exposure time) the 
number of counts available allowed us to divide this clus- 
ter only in three radial bins. The temperature profile of 
A1995 is consistent to be flat within the errors with a 
temperature around 9 keV. The abundance profile seems 
to have a positive gradient in the outer bin, however the 
errors are large and this increase in Z is not very signif- 
icant. 

3. SELF-SIMILARITY OF RADIAL PROFILES 

One of the main goals of this paper is to look for a 
(purely phenomenological) self-similarity in the radial 
profiles of temperature and metallicity, after they are 
scaled to the cluster virial radius rigo . A measure of rigo 
is thus crucial to test for such self-similarity in our clus- 
ter sample. This quantity can be approximated by the 
following relation: 

r 180 = 1.95/r 1 Mpc((kT)/l0 keV) 1/2 , (1) 

as calibrated from the non-radiative hydrodynamical 
simulations of clusters by Evrard et al. (1996). It is 
worth noticing that this relation is in agreement with 
the scaling relations observed (e.g. Ettori et al. 2004) in 
the X-rays, where the dependency on kT is consistent 
with Eq.(l) and only the absolute normalization may 
experience some variations. To compute the global 
temperature {kT) necessary to estimate rigo we extract 
spectra including emission going from 0.07rigo to 0.4rigo 
in each cluster. The central regions of each cluster are 
therefore excluded from the spectra in order to avoid 
contamination from a possible cool core. The values 
of (kT) and rigo have been evaluated iteratively until 
a convergence to a stable value of the temperature 
is obtained (AkT < 0.01 keV between two different 
iterations). From the fits we are able to determine also 
a global metallicity (Z) in each cluster. In Table [2] 
we list the best- fit values for (kT) and (Z), and the 
value of rigo computed using the formula above. The 
minimum and maximum apertures used to extract the 
total spectrum are listed as well. 



3.1. The temperature profiles 

Figure [2] shows the normalized temperature profiles for 
all the CC clusters (panel b) compared with the NCC 
clusters (panel d) . This figure has been obtained by nor- 
malizing the temperatures in each cluster to its average 
temperature (kT) computed from the total cluster spec- 
trum excluding the central 0.07rigo. The error- weighted 
mean and the best-fit results after fitting with single 
power-laws Y oc r^ are presented in Table [3l 

Within O.lrigo, the temperature profiles in CC objects 
increases with a slope [i — 0.25. Moving outwards, be- 
tween O.lrigo and the outer radial limit of our spectral 
analysis at w 0.5rigo, these profiles behave as r~ 0A . Non- 
cooling-cores systems have, on average, a profile that is 
almost flat at r < O.lrigo and then decreases rapidly as 





Fig. 4. — Error- weighted mean temperature profile of the CC 
(filled circles) and NCC (open diamonds) clusters in our sample at 
intermediate redshift (top panel). (Bottom) Comparison between 
our results and the local measurements in DM02 (shaded regions). 
The la errors on the means are plotted as solid lines (dark gray 
region for the local estimates) while the scatter (rms) in each data 
bin is shown as a dotted line (light gray region for local estimates). 



r -0 - 3 . In the outskirts, the temperature profiles of CC 
and NCC clusters show a significant discrepancy between 
their slopes, being NCC more deviant from the isother- 
mal case. 

Our best fit functional for the CC sample is fully 
consistent with the best fit functional form found by 
Vikhlinin et al. (2005) in their sample of CC clusters, 
at r < 0.3rigo. The two functionals diverge significantly 
only above 0.3rigo, where our profile is flatter (and there- 
fore the value of kT j (kT) is higher) than Vikhlinin et al. 
(2005) profile. However, only a few of our data points are 
located beyond that radius, preventing us from any sta- 
tistically significant comparison between the two samples 
at r > 0.3rigo. 

Adopting the A05 abundances has not changed the 
best-fit values of kT at every radius (always fully con- 
sistent within the la errors) in the individual clusters, 
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TABLE 3 

Error-weighted means, with errors on the mean and r.m.s. quoted within round brackets, and best-fit parameters of the 
SINGLE POWER-LAWS Y = Yq.1 (x/O.l)^ , WITH x = r/r^so- (kT) and (Z) are measured in the radial range 0.07 — 0.4rigo- 



kT/{kT) 


All 


all r 
0.84 ±0.04(0.28) 
Y i = 1.00 ±0.01 

A* = 0.15 ±0.01 
x 2 /dof=1223.4/78 


r < O.lriso 
0.79 ± 0.03(0.35) 
Y i = 1.20 ±0.01 
At = 0.26 ± 0.01 
x 2 /dof=720.0/43 


r > O.lriso 
1.02 ±0.06(0.10) 
Y i = 1.10 ±0.02 
At = -0.16 ±0.03 
x 2 /dof=59.1/33 


CC 


0.80 ± 0.04(0.24) 
y .i = 1.00 ±0.01 
fj, = 0.17 ± 0.01 
x 2 /dof=744.3/52 


0.76 ±0.03(0.26) 
Y .i = 1.16 ±0.02 
At = 0.25 ± 0.01 
x 2 /dof=496.3/30 


1.02 ± 0.06(0.08) 
Y .i = 1.07 ±0.02 
At = -0.10 ± 0.03 
x 2 /dof=26.7/20 


NCC 


1.13 ±0.07(0.20) 
y .i = 1.12 ±0.01 
At = -0.15 ±0.01 
x 2 /dof=76.6/24 


1.27 ±0.08(0.14) 
y .i = 1.25 ±0.04 
fi = -0.02 ± 0.03 
x 2 /dof=28.8/ll 

Z/(Z) 


1.02 ±0.07(0.13) 
Y .i = 1.23 ±0.04 
fi = -0.32 ± 0.05 
x 2 /dof=17.1/ll 


All 


all r 
1.31 ±0.29(0.54) 
Y i = 1.13 ±0.04 
fi = -0.23 ±0.02 
x 2 /dof=149.2/78 


r < O.lriso 
1.41 ±0.26(0.60) 
Y i = 1.03 ±0.06 
fi = -0.28 ±0.03 
x 2 /dof=120.7/43 


r > O.lriso 
1.11 ±0.33(0.32) 
Y i = 1.35 ±0.10 
^ = -0.44 ±0.15 
x 2 /dof=20.0/33 


CC 


1.40 ±0.28(0.58) 
Y i = 1.17 ±0.05 
At = -0.24 ±0.03 
x 2 /dof=80.2/52 


1.53 ±0.26(0.62) 
Yoi = 1.11 ±0.06 
At = -0.27 ±0.03 
x 2 /dof=64.9/30 


1.12 ±0.32(0.33) 
Y i = 1.37 ±0.11 
fi = -0.52 ±0.18 
x 2 /dof=l 1.0/20 


NCC 


1.06 ± 0.31(0.44) 
Y i = 1.06 ±0.07 
At = 0.00 ± 0.07 
X 2 /dof=41.7/24 


1.05 ±0.28(0.56) 
Y i = 1.08 ±0.15 
At = 0.03 ± 0.14 
x 2 /dof=31.7/ll 


1.08 ±0.36(0.31) 
Y i = 1.27 ±0.20 
fi = -0.29 ±0.27 
x 2 /dof=8.5/ll 



therefore the best fit functionals representing both the 
CC and the NCC sample have not varied. 

A clearer picture can be seen also if we compute an 
error-weighted average of the kT profiles in several bins 
of width 0.05 in r/r 18 o- The contribution to the single 
bin is provided from the measurements (and relative er- 
ror) that fall into that bin, weighted in proportion to the 
percentage of the spatial coverage of the bin. The error- 
weighted mean kT j (kT) profile is plotted in Figure H] 
and compared with the local estimates from DM02. CC 
and NCC objects show well defined opposite gradient 
in the inner radial and more similar behaviour moving 
outwards. A good agreement is also observed with the 
local profiles, apart from two significant deviations: (i) 
our CC mean profile appears flatter at r > 0.2riso, with 
an error-weighted value of 0.97 ± 0.06 (r.m.s. 0.07), to 
be compared with the local value of 0.85 ± 0.11 (r.m.s. 
0.16), (ii) our NCC profile is steeper within O.lrigo, with 
a mean value of 1.27 ± 0.08 (r.m.s. 0.14) with respect to 
the local value of 1.04 ± 0.05 (r.m.s. 0.07). 

3.2. The metal abundance profiles 

The metallicity profiles are plotted against the radius 
normalized to ri§o for the CC and NCC sample in Fig- 
ure [3)3 and[3Ji, respectively. A different behaviour in the 
very central regions between the two samples is quite 
clear. To characterize this behaviour, we have fitted the 
normalized profiles Z/ (Z) with respect to r/r igo with sin- 
gle power-laws Far' 1 over different radial ranges. While 



the NCC clusters presents a flat profile within ~ O.lrigo, 
a sharper negative gradient is observed in the CC cluster 
sample (^ = —0.24; see Table [3j. Also at r > O.lrigo, 
where a model with a single power-law well reproduce 
the data (reduced % 2 less than 1), the CC clusters show 
hints of a steeper profile (a^ = —0.52 ± 0.18 in the CC 
sample, [i = -0.29 ± 0.27 in the NCC sample). 

We compute an error-weighted average Z profile as 
done for the temperature profile and compare it with the 
local measurements in DM01, after scaling the radii to 
rigo (see Figure [5]). The two profiles are in agreement at 
r > O.lrigo, with both subsamples showing an evidence 
for a negative gradient in metallicity significant at least 
at 2er. A different behaviour between the two subsam- 
ple is observable in the central bin, where the value in 
the CC sample is - 20% - 30% higher than in the NCC 
sample. It is worth noticing that all the clusters in our 
sample have kT > 6 keV and this trend may be differ- 
ent in lower temperature clusters. DM01 in their analy- 
sis observed a clear gradient in the metallicity profiles of 
their CC clusters, while the profiles of their NCC clusters 
were almost constant. Moreover the average metallicity 
observed in the CC clusters was systematically higher 
than in the NCC sample, at least within ~ 0.3r/rigo 
from the center. In our analysis, we do not find a clear 
difference between the CC and NCC abundance profiles 
as in DM01. Except for the inner radial bin where the 
metallicity in CC objects can be higher by 50 per cent 
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Fig. 5. — Error-weighted mean metal abundance profile of the CC 
(filled circles) and NCC (open diamonds) clusters in our sample at 
intermediate redshift (top panel). (Bottom) Comparison between 
our results and the local measurements in DM01 (shaded regions). 
The 1<t errors on the means are plotted as solid lines (dark gray 
region for the local estimates) while the scatter (rms) in each data 
bin is shown as a dotted line (light gray region for local estimates). 



than in NCC ones, the CC and NCC profiles are very 
similar and consistent within the errors. 

To compare our mean values at intermediate redshift 
with the results obtained locally from DM01, we esti- 
mate an error-weighted average Z profile for both CC 
and NCC objects in local and intermediate z samples 
(Figure [5|) . While the slope of the profiles is generally 
in agreement with DM01 for the both the CC and NCC 
clusters, the value of Zj (Z) is systematically higher in 
our sample with respect to the DM01 sample, with dif- 
ferences up to 50% within O.lrigo of CC systems. This 
might be due to the different method used to compute 
(Z), being in DM01 work estimated as fit with a con- 
stant to the radial metallicity profile. However, apart 
from the inner radial region, the discrepancy between 
the local profiles and the ones at intermediate redshits 



are within la. Using the same method to determine (Z) 
on our data, we find that the values of Z/(Z) arc fully 
consistent with DM01. 

3.3. Comparison with the new compilations of solar 

values 

All the analysis described in the current section has 
been performed adopting the AG89 compilation of pho- 
tosphcric abundances. This choice has been due mainly 
by the necessity to have a direct comparison with pre- 
vious works in the literature (e.g. DM01). However, 
as explained above in §1, the abundance values listed in 
AG89 have been recently superseeded by the new pho- 
tosphcric values by Grevesse & Sauval (1998) and A05, 
who introduced a 0.676 and 0.60 times lower Iron solar 
abundance, respectively. Therefore, we have also per- 
formed the fits using solar abundances by A05 to check 
whether adopting the "old" AG89 values might have in- 
troduced any bias in our analysis. 

The shape of the Z profile for both the CC sample and 
the NCC sample resembles very closely that observed in 
Figure [5] for the AG89 values of Z. However to better 
quantify this comparison we fitted also these new values 
with a power-law functional. We measure 

Z CC /(Z) = 0.67 ± 0.06 s-o.25io.03. (2) 



Zncc/(Z) = 1.00™ ar ™- 07 , (3) 

with x = r/risQ. If we compare the result of the fit with 
the last two rows in the first column in Table [3l it is 
quite clear that the slope of the power-law in both the 
CC and NCC sample is consistent with what we obtained 
using the AG89 values (well within the la uncertainties) 
and the difference is only in the normalization. As ex- 
pected, this result might indicate that our metallicities 
are mostly driven by Iron and that the contribution of 
the a elements to the determination of the abundances 
is negligible. Indeed, all our clusters have a temperature 
larger than ~ 6 keV, therefore the abundance measures 
are dominated by the Fe-Ka line. As a further test to 
this hypothesis we tried to fit the Fe abundance indepen- 
dently from the a elements abundances. To this aim we 
used a vmekal model where the abundances of O, Mg, Si 
and S were tied-up to the same value and fitted as a sin- 
gle free parameter (in order to reduce the number of free 
parameters), while the other elements (apart from Fe) 
were frozen at solar. While the Fe abundance remained 
always consistent with the value of Z measured consid- 
ering the metallicity as a single parameter, in almost all 
the spectra we could not find any statistically significant 
detection of a contribution from a elements. Their abun- 
dances has been measured as upper limits in most cases 
(generally Z a < 0.3), and as very low values in the rest 
of the spectra (generally Z a ~ 0.1 — 0.2, often consistent 
with Z a ~ at la). This is true even in the inner part 
of the clusters in the sample, where the signal-to-noise of 
the spectra is higher and in principle it may be easier to 
detect the presence of elements other than Iron. 
These results suggest that the measure of Z in our clus- 
ter sample consists mainly in a measure of Iron metal- 
licity. Therefore, adopting the AG89 solar abundances 
instead of the A05 (which have different Zp e /Z a ratios) 
produces only a difference in the value of the relative Z 
measured. This does not introduce any bias in the anal- 
ysis of the radial profiles, since the absolute value of Z 
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Fig. 6. — Distribution of the exponential fi in the relation 
kT/(kT) = A (r/riso)' 1 a s a function of T c /r age . 



does not change and only the reference value assumed 
for the solar metallicity experience a variation. 

3.4. Gradients and cooling times 

We investigate the correlation of the central slopes of 
the temperature and metallicity profiles with the cen- 
tral cooling time in each cluster. To this aim, the val- 
ues of kT(r)/{kT), Z(r)/(Z) and r/ri 80 (normalized to 
the average temperature, the average metallicity and the 
virial radius measured in each cluster, respectively, as 
described in Section 3) with r < O.lriso have been con- 
sidered to characterize the cooling cores. We find that 
the most robust correlation is present between the slope 
p of the temperature profiles, 

kT/(kT) = A(r/r 180 r, (4) 

and Tc/rage, with a Spearman's p rank correlation value 
of -0.87 that corresponds to a significance of the non- 
correlation case of P = 2 x 10~ 4 . In the kT — Z and 
Z — r relations, the values of the Spearman's p are 0.44 
and 0.54, corresponding to a significance of 0.15 and 0.07, 
respectively. The exponent p correlates with r c /r age , be- 
ing higher at lower values of T c /r age , with all the CC 
clusters having p > and r c /r age < 0.6 (Figured]). 

4. CONCLUSIONS 

In the present work we analyzed a sample of 12 
galaxy clusters present in the Chandra archive with 
at least ~ 20, 000 net ACIS counts and kT > 6 keV. 
These clusters were chosen in the 0.1-0.3 redshift range, 
regardless of their shape. We computed the cooling time 
of the clusters, subdividing the sample in 7 cool core 
clusters and 5 non-cool core clusters. This subdivision 
allowed us to compare the two categories in a systematic 
fashion, following the approach of DM01. We performed 
a spectral analysis in radial bins of each cluster in the 
sample requiring each bin to have ~ 7000 — 8000 counts, 
fitting the spectra with a thermal model with Galactic 
absorption. This allowed us to derive temperature and 
metallicity profiles for each cluster. The virial radius 



ri8o was computed in order to renormalize the radii 
to physically meaningful quantities and investigate for 
self-similarities in the radial profiles. To this aim the 
global temperature {kT) and metallicity (Z) in each 
cluster were measured as well. The main results coming 
from our work can be summarized as follows. 

• The temperature profiles in the inner O.lriso have, 
on average, a positive gradient, kT{r) oc r M with 
p ~ 0.25 in CC systems, whereas it is almost flat 
in NCC systems. The outer regions are well fitted 
with a single power-law with slopes significantly 
different, being steeper {p = -0.32 ±0.05) in NCC 
objects. The general trend of our CC sample is 
fully consistent with Vikhlinin et al. (2005) at r < 
0.3ri8o. The low number statistics above 0.3riso 
prevents us from any statistically significant com- 
parison between the two samples at r > 0.3riso- 

• The metallicity profiles in the inner regions is al- 
most constant in NCC clusters around the value 
measured excluding counts from r < 0.07riso- In 
the CC sample, a steep negative gradient is ob- 
served (p — —0.27 ± 0.03) in the central regions. 
At r > O.lriso, a power-law reproduces well the 
distribution of the spectral measurements, with 
a slope that is marginally steeper in CC clus- 
ters (p = -0.52 ± 0.18) than in NCC clusters 
(p = -0.29 ±0.27). 

• Comparing our averaged metallicity profiles with 
the ones in DM01, we found that our values of 
Zj (Z) are systematically higher, with differences 
up to 50% within O.lriso of CC systems. This may 
be explained by the different method adopted in 
DM01 to estimate (Z), as best- fit with a constant 
over the entire metallicity profile, without any ex- 
clusion of the central core. 

• Using the solar abundances from Asplund (2005, 
A05) gives consistent results with what we ob- 
tain using the values by Anders & Grevesse (1989, 
AG89), with a discrepancy only in the normaliza- 
tion (as expected, ~ 60 — 70% higher) but not in 
the slope of the Z radial profiles. Together with the 
fact that, in most cases, we were able to measure 
the a elements only as upper limits, this indicates 
that our metallicities are mostly driven by Iron and 
that adopting the AG89 solar abundances instead 
of the A05 results in a difference only in the abso- 
lute values of the Z measured but does not intro- 
duce any bias in the radial profile analysis. 

• Fitting a power-law shape to the temperature pro- 
files, kT/\kT) = A (r/ri 80 ) M , we found that p cor- 
relates strongly with the cluster cooling times, be- 
ing higher at low values of T c /T age , with all the CC 
clusters having a > and r c /r age < 0.6. As ex- 
pected, strong correlation is also observed between 
the inner slope of the metallicity profile and cluster 
cooling time. 

In general, our results further demonstrate the in- 
valuable role played by X-ray archival studies of the 
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chemo- and thermo-dynamical properties of galaxy clus- 
ters. Analyses based on the Chandra archive, like that 
presented here (see also Vikhlinin et al. 2005; Balestra 
et al. 2007; Maughan et al. 2007), in combination with 
analogous studies from the XMM-Newton archive, will 
constitute an important heritage from the present gen- 
eration of X-ray satellites for years to come. Nowadays, 
available data on the evolution of the chemical enrich- 
ment of the ICM provide important constraints on mod- 
els aimed at explaining the past history of star formation 
and the dynamical processes taking place during the cos- 
mological build up of galaxy clusters. However the study 
of the thermo-dynamical properties of the cooling cores 
and the evolution of the abundance distributions in clus- 
ters with the redshift are just a part of what can be cur- 
rently done exploiting in full the existing Chandra and 
XMM-Newton archives. Archival works like ours have 



the potential to shed new light on the properties of the 
stellar populations responsible for the ICM enrichment, 
and on the mechanisms which lead to the generation of 
the cool cores and determine the transport and diffusion 
of heavy elements from star forming regions. 
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